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ABSTRACT 

Observations of HDO are an important complement for studies of water, because they give strong constraints on the formation processes - grain 
surfaces versus energetic process in the gas phase, e.g. in shocks. The HIFI observations of multiple transitions of HDO in Sgr B2(M) presented 
here allow the determination of the HDO abundance throughout the envelope, which has not been possible before with ground-based observations 
only. The abundance structure has been modeled with the spherical Monte Carlo radiative transfer code RATRAN, which also takes radiative 
pumping by continuum emission from dust into account. The modeling reveals that the abundance of HDO rises steeply with temperature from 
a low abundance ( 2.5 x 10~ H ) in the outer envelope at temperatures below 100 K through a medium abundance ( 1.5 x 10~ 9 ) in the inner 
envelope/outer core, at temperatures between 100 and 200 K, and finally a high abundance ( 3.5 x 10~ 9 ) at temperatures above 200 K in the hot 
core. 

Key words. ISM: abundances — ISM: molecules 



1. Introduction 

Water is known to be a fundamental ingredient of the interstel- 
lar medium. It is a major coolant of star-forming clouds, but, 
apart from some high-lying maser transitions, is unobservable 
from the ground because of the atmospheric absorption. The 
satell ites SWAS and Odin dNeufeld et alll2003tlSandqvist et all 
2006) have observed the ground-state orth o-water line, and 
ISO has observed highly excited transitions (Goicoechea et al., 
200-4 ICernicharo et all l2006t [Polehampton et al., 2007)), al- 
though with low spectral resolution. Water is also an important 
reservoir of oxygen and therefore a key ingredient in the chem- 
istry of oxygen-bearing molecules. Investigations of the water 
abundance using many lines covering a wide range of excita- 
tion levels with high spectral resolution have been among the 
main reaso ns for building the HIFI instrument on board Herschel 
' dde Graauw et aUl2010UPiibratt et aUl2010h . 

Water can be formed in the gas phase through an ion- 
molecule channel that is initiated by a reaction of H + or W\ with 
atomic oxygen. Alternatively, water can also be formed through 
highly endothermic reactions that involve atomic oxygen, which 
reacts with molecular hydrogen to form OH, which in turn can 
react with molecular hydrogen to form H2O. The first process 
can work in rather diffuse clouds, where the ionization fraction 
is high, while the second one is efficient only in shocks, because 
only there the temperatures are high enough to overcome the 



* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



high endothermicity. In dense gas, which is neither highly ion- 
ized nor very hot, these channels do not work efficiently. There it 
is assumed that water is formed on grain surfaces, where oxygen 
can react with atomic hydrogen. These icy grain mantles can be 
studied through obse rvations of t he wa ter stretching bands in ice 
in the infrared (e.g., iGibb et all |2000), by which the existence 
of water ice is a well established observational fact. It is then 
assumed that this ice mantle is returned to the gas phase either 
through thermal heating (where water ice mantles are destroyed 
at about 100 K), or through shock sputtering. The net result is 
that gas-phase water, irrespective of its formation channel, is of- 
ten observed in regions where the energetics are such that it is 
not a priori clear if the water has a gas phase or grain surface 
origin. 

Deuterium fractionation, i.e. the enrichment of the deuter- 
ated counterpart of a molecular species much above the cos- 
mic [D]/[H] ratio, is a process that is caused by the zero-point 
energy difference of the vibrational potential. The deuterated 
species have a lower zero-point energy because of their larger 
reduced mass. For H2O and HDO, the difference in energy is 3.4 
kJ/mole or 409 K. That means that deuterium fractionation for 
H2O takes place mainly at low temperatures. In other words, the 
[HDO]/[H20] ratio should be high if the water in question has 
been formed at a low temperature, i.e. on cold grains, but low if 
it is a product of either PDR or shock chemistry. Observing the 
[HDO]/[H20] ratio is therefore an important clue to the forma- 
tion mechanisms of H2O 

This was the goal of the study by IComito et al.l |2003, 
hereafter C03), who investigated the ground-state HDO line at 
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893 GHz observed with the Caltech Submillimeter Observatory 
towards Sgr B2(M) and (N), as well as some published high- 
excitation J1DO linesobserved with the IRAM 30-m tele- 
scope dJacq et all [1990). The HDO observations were com- 
pared with those of the o-Hl, 8 ground-state line from SWAS 
dNeufeld et all 120031) . and of hig her-energy p-Hj 8 Q trans itions 
observed with the IRAM 30-m bv lGensheimer et al.1 (Il996t) . 

The HDO and Fli 8 ground-state lines, observed in absorp- 
tion, are formed in the outer envelope, while the high-energy 
transitions trace the central hot core. Up to now, no observations 
tracing the inner Sgr B2(M) envelope/outer core were available. 
IComito et al. (2003) concluded that in the outer layer, H2O was 
predominantly shock-produced, while the [HDO]/[H20] ratio in 
the hot core showed, as expected, a grain surface origin. Recent 
HIFI observations in the framework of the Guaranteeed Time 
Key Program Herschel Observations of Extra-Ordinary Sources 
(HEXOS) make it now possible to extend the study of the HDO 
distribution to the full envelope of Sgr B2(M). 

Sgr B2(M) and its neighbor Sgr B2(N) are some of the most 
massive star-formation sites in our Galaxy. Understanding their 
history and structure will help understanding the formation of 
massive Galactic center clusters, such as Arches, and will also 
have ramifications for even more extreme star-formation pro- 
cesses in starburst galaxies, and even in the early universe. 

HDO level diagram 



Table 1. Observed HDO lines used in the modeling. 




Fig. 1. Level diagram of the HDO lines in the present dataset. 
The frequency of each transition is indicated in GHz. Transitions 
are color-coded as follows: in solid red, HIFI observations 
in band la and lb (this work); in dashed orange, transi- 
tions falling within the overall HIFI range (to be observed 
within the HEXOS Key Program); in dark blue, submillimeter- 
wavelength observations (C03; M. Gerin, priv. c omm.); in 
light b lue, millimet er- wavelen gth observations from Jaca et all 
(1990), Nummelin et al.l d 19981). and Belloche et al., in prepa- 
ration (cf. Belloche et al., 2008); in black, transitions that fall 
outside of the HIFI range. 



2. Observations 

Full spectral scans of HIFI bands la and lb towards Sgr B2(M) 
(«72000 = 17 /l 47 m 20.35 s and <5, 2 ooo = -28°23'03.0") have been 
carried out on March 20101 and 2, providing coverage of the 
frequency range 479 through 637 GHz. The HIFI Spectral Scans 
are carried out in Dual Beam Switch (DBS) mode, where the 
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Rows 1 through 5: summary of all the HDO transitions observed in the 
HIFI bands la and lb towards Sgr B2(M) (this work). Rows 6 through 
9: list of previously observed HDO features that were used to constrain 
the results of our model (see Sect. F3J. References: a), this work; b), 
Belloc he et al., in preparation; c). lJacq et aljdl990i): < j).lNummelin et all 
dl998h : e), M. Gerin, priv. comm.: f) JComito et all J2003I) . 



DBS reference beams lie approximately 3' east-west or 6' apart. 
The wide band spectrometer (WBS) is used as a back-end, pro- 
viding a spectral resolution of 1.1 MHz over a 4-GHz-wide 
Intermediate Frequency (IF) band. A HIFI Spectral Scan con- 
sists of a number of double-sideband (DSB) spectra, tuned at 
different Local Oscillator (LO) frequencies, where the spacing 
between LO settings is d etermined by the "redundancy" selected 
by the observer (Comito & Schilkg, l2002l) . The molecular spec- 
trum of Sgr B2(M) in band la and lb has been scanned with 
a redundancy of 4 and 8, respectively, which means that every 
frequency has been observed 4 and 8 times respectively in each 
sideband. Multiple observations of the same frequency at dif- 
ferent LO tunings are necessary to separate the lower-sideband 
(LSB) from the upper-sideband (USB) spectra. 

The data were calibrated thro ugh the sta ndard pipeline dis- 
tributed with version 2.9 of HIPE dOti,l2010t) . and subsequently 
exported to CLASSQ using the HiClass task within HIPE. 
Deconvolution of the DSB data into single-sideband (SSB) for- 
mat was performed with CLASS. All the HIFI data presented 
here, spectral features and continuum emission, are deconvolved 
SSB spectra. Both horizontal (H) and vertical (V) polariza- 
tions were obtained. A discrepancy in flux of up to 20% is ob- 
served between the two. Because the H data in band la and 
lb show fewer instrumental artifacts than the V data, we will 
here show H-polarization spectra only. The intensity scale is 
main-beam temperature, and is a results of applying a beam ef- 
ficiency correction of .69 for band la, and 0.68 for band lb 
(Roelfsem a et alluOlOl) to the deconvolved data . 

Table Q] and the level diagram in Fig. [TJ summarize the 
HDO transitions that were observed so far towards Sgr B2(M). 
Together with those observed with HIFI (top five lines in TableQ] 
red transitions in Fig. [TJ, our dataset includes mm-wavelength 
transitions from the IRAM-30m telescope and from SEST, and 
submm-wavelength transitions acquired with the CSO (refer- 
ences are indicated in Tab. [TJ. The frequencies were obtained 
from the JPL Molecular Spectroscopy Catalog dPickett et all 
1998). 



1 Continuum and Line Analysis Single-dish Software, distributed with 
the GILDAS software, see http://www.iram.fr/IRAMFR/GILDAS 
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3. Results 

The spherical Monte Carlo radiat ive transfer code RATRAN 
(Hogerheijde & van der Tak, 2000) is employed to compute the 
intensity of the molecular lines and of the dust co ntinuum. The 
sourc e model (for a more detailed description see iRolffs et all 
2010) corresponds to a centrally heated source with a density 
distribution that follows a radial power law with an index of 
1.5. The best match between the observed and modeled contin- 
uum is achieved through the basic MRN grain size distribution 
(Ossenkopf & Henning, 1994, and references therein), that is, no 
formation of ice mantles and no coagulation onto the grains is 
included. The assumed source luminosity is 6.3 x 10 6 Lq. At 
the inner radius, 1365 AU, the density is 1.5 x 10 8 cm" 3 , and 
the temperature is 1400 K. The temperature is computed in an 
approximate way from central heating, using the diffusion equa- 
tion in the inner part and the balance between heating and cool- 
ing in the outer part, with an effective temperature of the dust 
photosphere of 54.5 K. The temperature has a steeper gradient 
in the inner part (200 K at 0.07 pc; 100 K at 0.14 pc) and de- 
creases to 15 K at the radius of 12.3 pc, where the density is 
2.7 x 10 3 cm -3 . From this point on, the temperature and den- 
sity remain constan t, until t he outer envelope radius of 22.5 pc is 
reached (iLis & Goldsmith!. 1 19891) . This model yields a total H 2 
column density of 7 x 1 24 cm -2 . The HDO collisional rates have 
been provi ded by LAMBDA (L eiden Atomic and Colli sional 
Database, ISchoier et all 120051) . based on the work of iGreenl 
<fl989h . 

Figure [2] shows in grey the three HDO emission features de- 
tected with HIFI, and also the CSO absorption feature at the bot- 
tom. Our best-fit model is overlaid in black. The observed HDO 
features display a half-power width (HPW) of 14.5 km-s" 1 , and 
are centered at \lsr - 63 km-s -1 , with the exception of the 894- 
GHz absorption feature, which is centered at Vlsr = 64.7 km-s -1 
(see Fig. |2|. This discrepancy is within the uncertainty on the 
velocity scale that affects HIPE 2. In this work we will use the 
same velocity parameters (Ylsr = 63 km-s -1 , Av= 14.5 km-s -1 ) 
for all features in the model. 

It has already been demonstrated (C03 and references 
therein) that a constant HDO abundance across the whole en- 
velope of Sgr B2(M) cannot reproduce the observed intensities 
of features spanning a relatively wide range of excitation, which 
means that they arise from different regions within the cloud. It 
is reasonable to expect a drop in abundance when the gas and 
dust temperature fall below 100 K, causing the water content of 
the gas to freeze out onto the dust grains. 

With the current model, a simultaneous fit of all the HIFI and 
CSO transitions (see Fig. Q]) is obtained assuming that [HDO] 
= 1.5 X 10 -9 when T gfls > 100 K, and [HDO] = 1.3 x 10 -11 
when T gas < 100 K. However, this two-step abundance model 
underestimates the intensity of the mm-wavelength features. 

Figure [3] shows the normalized populations of the HDO 
levels of the relevant features in our dataset as calculated by 
RATRAN. The bulk of the emission for the J=3 U - 2 2 ,i 226- 
GHz transition is produced in the inner hot-core region, where 
T> 200 K. The other features though arise predominantly in the 
inner-envelope/outer-core region, with temperatures close to or 
lower than 100 K. In order to reproduce the intensities of the 
3 1,2 - 22,i transition as well, we need to introduce a rise in HDO 
abundance to about 3.5 x 10 -9 , where T gas > 200 K. 

Our model very precisely reproduces 8 out of the 10 transi- 
tions in our sample (see Tab.|2}. On the other hand, the 1 i,o - 1 1,1 
and the 2\\ —1\^ transitions, at 80.6 and 241.6 GHz respectively, 
are underestimated by about a factor 3. Although we cannot 



Table 2. Observed and modeled intensities for all the features 
in our dataset (see TableQ}. 
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completely rule out contamination by other species, especially 
at 3-mm wavelength, where the contribution of unidentified fea- 
tures to the spectrum is about 50% (Belloche et al., 2008), we 
do not think it realistic that these two features are blended with 
other, weaker ones to the point of appearing three times stronger 
than they are. The model's failure to reproduce these lines may 
be due to the assumption of spherical symmetry, which does not 
capture all the details of the actual source structure. 



4. Discussion 

With the Herschel lines discussed in this paper, we now pos- 
sess information about the HDO abundance throughout the en- 
velope (Fig. 0. The physical structure model we employed is 
based on the struct ure used bylComito et alj d2003l) . refined by 
the study of IRolffs etall d2010l) . We note that we can fit 8 out 
of 10 lines using this model and a very simple abundance distri- 
bution: a low abundance (1.3 x 10 -11 ) below 100 K, a medium 
abundance (1.5 x 10 -9 ) at temperatures between 100 and 200 K, 
and a high abundance (3.5 x 10~ 9 ) at temperatures above 200 K. 
This is in contrast to the simpler two-step abundance structure 
in C03, where no constraints on the HDO abundance from the 
inner envelope/outer core were available. The exact boundaries 
of these abundance changes are not well determined, and while 
evaporation of ice produces a real discontinuity, other processes 
do not. As mentioned in Sect. Q] no process but ice evapora- 
tion can produce HDO in noticeable amounts. In the hot gas 
phase, this molecule will be formed relative to water with the 
same abundance of deuterium relative to h ydrogen, 1 .7 x 10~ 6 in 
the Galactic center (Lub owich et al.ll2000h which, even if all the 
oxygen were tied up in water, would correspond to an abundance 
of only a few times 10 -10 . This is almost two orders of magni- 
tude lower than what is observed in the hot gas. A possible ex- 
planation for the abundance gradient in gas at T > 100 K could 
be that the ice evaporation is not complete at 100 K, and that 
the gas-phase HDO abundance gradually increases until a tem- 
perature of about 200 K is reached; or that in the inner, denser 
regions, simply more water (and therefore HDO) has formed on 
the ice mantles. 

To make use of the full information contained in this study, as 
described in Sect.[T] we need to complement it with the study of 
H2O and its isotopologues, and then determine the [HDO]/[H20] 
ratio throughout the envelope. The water modeling will not be as 
straightforward, because the higher opacities will make the H2O 
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Fig. 2. Grey: the top three panels show the HDO emission fea- 
tures detected with HIFI, the bottom panel shows the ground- 
state absorption feature detected with the CSO. Our best-fit 
model, described in § [3] is overlaid in black. Note that the ab- 
sorption line is about 2 km-s -1 off from the center velocity of 
the emission lines, which is within the uncertainty of the HIFI 
calibration in HIPE 2.9. 



transitions much more sensitive to deviations from the spherical 
symmetry. Moreover, water will be produced in outflows (within 
the core) and shocks (in the outer envelope). However, the com- 
parison will provide necessary clues to the study of the develop- 
ment of the Sgr B2(M) envelope. 
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